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ABSTRACT 


This report presents tables of elementary three-dimensional absorbed- 
dose distributions in a water phantom irradiated by monoenergetic, 
point-monodirectional electron beams. Such distributions have been 
Obtained by the Monte Carlo method for 14 beam energies from 1 MeV 

to 60 MeV. The tabulated results can be applied to the determination 
of absorbed-dose distributions from parallel beams of arbitrary finite 
cross section. The beam of interest is treated as a superposition of 
point-monodirectional beams, and the absorbed-dose distribution is 
obtained as a corresponding superposition of elementary absorbed-dose 
distributions. By way of example, the tabulated data are used to 
obtain (1) depth-dose curves and practical ranges in broad-beam 
geometry, and (2) central-axis depth-dose curves, radial distributions 
of absorbed dose, and isodose patterns for beams with finite circular 
cross section. 


Key Words: Dosimetry; electrons; Monte Carlo; point-monodirectional 
beams; Superposition; treatment planning. 


This report is based in part on a paper presented at the Electron 
Dosimetry and Arc Therapy Symposium, Wisconsin Clinical Cancer Center, 
University of Wisconsin, Madison, Wisconsin, September 10-11, 1981. 


1. INTRODUCTION 


The determination of absorbed-dose distributions from electron beams, with the 
accuracy required by treatment planning, is a complex problem whose solution is 
facilitated when experimental dosimetry is supported by calculations. Flexible, 
approximate, analytical methods are required as well as detailed transport 
calculations. Among the latter, the Monte Carlo method (track simulation) is 
probably the most expeditious if not the most economical approach, because it 
requires the fewest assumptions and approximations and allows the fullest 
utilization of the available cross sections for the interactions of electrons 
with matter. The Monte Carlo calculations fall into two categories. The first 
category consists of special-purpose calculations that simulate as closely as 
feasible the beam and medium configurations encountered in practice. The second 
category consists of calculations for schematized and somewhat oversimplified 
configurations, and is useful for general surveys and parametric studies; it 
also can provide benchmark results for testing experimental dosimetry methods, 
and may be useful for checking the validity of assumptions made in approximate 
analytical treatments. 


The present paper describes Monte Carlo calculations that fall into the second 
category. Attention is focussed on the application of the superposition method 
[1-4]* to the calculation of absorbed-dose distributions in a homogeneous water 
phantom. According to this method, the electron beam of interest is treated as 
a superposition of point-monodirectional beams, and the spatial distribution of 
absorbed dose is obtained by a corresponding summation over the elementary 
absorbed-dose distributions from point-monodirectional beams. 


Such elementary distributions have been obtained for monoenergetic point- 
monodirectional beams with energies from 1 MeV to 60 MeV. Some of the character- 
istics of these elementary distributions are described here, and examples are 
given of the application of these results to other beam geometries. 


2. ELEMENTARY ABSORBED-DOSE DISTRIBUTIONS 


2.1. Definitions. The phantom is assumed to be a homogeneous semi-infinite 
medium that occupies the region z2>0. A point-monodirectional (i.e., a narrow 
pencil) beam of electrons with kinetic energy To, traveling along the z-axis, 
is assumed to be incident perpendicularly” onto the surface of the phantom. The 
Monte Carlo method has been used to obtain the elementary distribution function 
F(p,z), defined such that F(p,z)dodz represents the average fraction of the 
incident beam energy that is deposited at depths between z and z+ dz and at 
radial distances from the z-axis between po and p+dpo. The absorbed dose at 


a point (p,z) is given by (To/d) F(e,z)/2nm9, where d is the density of the 
medium. 


‘Figures in brackets indicate literature references at the end of this paper. 


*The assumption of perpendicular incidence is appropriate for electron beams that 
have been broadened by scanning, but does not take into account the angular 
dispersion of a beam broadened by a scattering foil. In principle it is possible 
to consider elementary point-monodirectional beams incident at various angles, 


and to treat the angular dispersion by an appropriate summation over such 
elementary beams. 


2.2. Method. The Monte Carlo model, which has been described earlier [5-7], 
takes into account the transport of energy by the primary electrons from the 
beam, and by all generations of secondary electrons and photons. The transport 
of electrons is simulated by dividing their tracks into many small segments, and 
sampling the multiple-scattering angular distribution in each successive segment 
from the Goudsmit-Saunderson distribution and the energy loss from the Landau 
distribution. The histories of bremsstrahlung photons emitted by the electrons 
are followed by sampling individual Compton-scattering, pair-production and 
photo-electric absorption events. The electrons set in motion in these inter- 
actions, and the secondary electrons resulting in ionization events, are treated 
in the same manner as the primary electrons. The energy deposition along the 
electron tracks is scored and classified to obtain the distribution F(p,z). 


2.3 Ranges. When graphing or tabulating elementary absorbed-dose distributions, 
Or other distributions derived from them, one finds it useful to express the 
distances o and z in units of the CSDA range? ro. When this is done, the 
distributions corresponding to different beam energies To "Scale" in the sense 
that the explicit dependence on Ty is minimized, which in turn facilitates 
interpolation with respect to To. 


In the treatment of energy loss by ionization and excitation in the Monte Carlo 
calculations whose results are reported here, the stopping properties of water 
were characterized by a mean excitation energy I = 65.1 eV, and the density 
effect (reduction of stopping power due to the polarization of the medium) was 
evaluated according to the prescription of Sternheimer and Peierls [8]. In 

recent work on stopping power we were led to the conclusion that improved 
accuracy would be obtained with a mean energy I = 75.0 eV, and with the density- 
effect correction derived by Ashley [9] from experimental optical data for water. 
The old ranges that correspond to the cross sections used in the Monte Carlo 
calculations are compared in Table 1 with the new ranges that would result from 
the use of a mean excitation energy I = 75.0 eV and Ashley's density-effect 
calculations. The new ranges are larger than the old ranges by 0.3% at 60 MeV, 
1.2% at 10 MeV and 1.7% at 1 MeV. Inasmuch as these differences are small, it 
does not appear necessary to repeat the Monte Carlo calculations with the 

improved cross-section input data. It seems adequate to assume that the absorbed- 
dose distributions, as functions of p/ro and Z/ro, remain unchanged, so that 
it is merely necessary to use the appropriate new range values when expressing 

the absorbed-dose distributions as functions of o and Zz. 


2.4. Elementary Absorbed-Dose Distributions from Point-Monodirectional Beams. 

A typical elementary distribution of absorbed energy from a 20-MeV point- 
monodirectional beam, is shown in Fig. 1. The information is equivalent to that 
contained in the distribution F(p,z) as defined above, but the representation 
is slightly modified. The histogram in the upper right-hand corner is a depth- 
dose curve in the form of a plot of the dimensionless quantity (ro/To)D vs. 
the scaled depth z/ro, where @ is the energy deposition per unit depth 


“The CSDA range ro, evaluated in the continuous-slowing-down approximation, 
represents a rectified path length over which the electron kinetic energy is 
reduced from To to zero. It is calculated as an integral with respect to 
energy over the reciprocal of the total stopping power. 


(obtained by integrating F(p,z) over all p). The other histograms show, for 
the six indicated (shaded) depth intervals, the radial distribution of deposited 
energy, as a function of the scaled distance p/rq from the z-axis. Each radial 
histogram bin represents the fraction of the energy deposited in the indicated 
ring around the z-axis. Each radial histogram is normalized to unit area. 


Figure 2 gives elementary absorbed-dose distributions in a different form, in 
terms of the normalized cumulative distribution 


0 CO 
2(0,z) - f F(o'sz)do'y f Flot sZido nee. (1) 
(6) O 


For point-monodirectional beams with energies of 5, 10, 20, and 40 MeV, contours 
are given as functions of the scaled depth Zz/ro along which © has the value of 
20, 40, 60, 80, 90, or 95%. For example, the curve of 99 vs. 2Z/Yro indicates 
the radial distance from the z-axis within which 90% of the energy is deposited 

at each depth. 


Figure 3 shows plots of %n(1-%) vs. (p/ro)?, at various depths in a medium 
irradiated with 10- and 20-MeV point-monodirectional beams. The quantity 1 - ® 
represents the fraction of the energy at a given depth that is deposited at 
radial distances greater than p/ro. If the radial distribution of deposited 
energy were Gaussian, an approximation made in certain approximate dosimetry 
calculations, then the curves in Fig. 3 would be straight lines. This is not the 
case, but a straight-line approximation might be adequate over the radial region 
around the z-axis in which ~ 3/4 of the energy is absorbed. 


2.5. Tables of Elementary Absorbed-Dose Distributions. Monte Carlo results for 
the distribution F(o,z) are given in Tables Al and A2 at the end of this paper for 
14 beam energies between 1 MeV and 60 MeV. At each energy, the spatial distribu- 
tion of deposited energy was obtained by analyzing the sampled electron-photon 
cascade from 10,000 incident electrons. Energy deposition was scored in annular 
volume elements. The depth and radial boundaries of these volume elements were 
expressed as fractions of the CSDA range ro, and the fractions were kept the 

same for each beam energy in order to take advantage of the scaling of the 
distributions. 


Table Al gives the distribution as a function of the depth only, 
Zi CO 
n 
H, ~ 100 f dz J ao Filo gz) 45s (2) 
ZA] 0 


where 2, =n ro/20; on -' |s2;...c0.. It Should be noted™ that \theyquamenty 
To H,/(Z_ - Zn-1) is an approximation to 100 9(z), where @(z) jis the energy 
deposition per iunit depthat iz = (2a eZ 2. 


Table A2 gives the normalized cumulative radial distribution 


Zz 
¥ (op) = 100 ee 5 (3) 


which is an approximation to the quantity 100 %(p,z) defined by Eq (1). 
3, SUPERPOSITION METHOD 


Let the incident electron beam be characterized by an intensity distribution 
Qo(Xo,¥o), such that Qo(x0,yo)dxodyo represents the number of point-monodirec- 
tional beams incident perpendicularly onto the surface of the phantom at lateral 
positions between Xo and Xo + dxo and between yo and yo + dyo. Let 
D(x,y,;Z) be the absorbed dose at the point (x,y,z), obtained by summing over 
the absorbed-dose distributions from the point-monodirectional beams. To 
Facilitate the summation, we replace the coordinates Xo, yo by cylindrical 
coordinates p,y, such that 


Qo(Xo02Yo)dxXodyo = Qo(x + pcosy,y + psiny)pdpdy = Q(p.y)pdody . (4) 


The absorbed dose can then be written as 


[ee] 


2 
p(xsysz) = 22 f odo fo “dy (ose) Fhes2) (5) 


0 O 


It is useful to introduce the beam-geometry function 


Zill 
G(x.y.0) = [ de Qo.) , (6) 
0) 


in terms of which the expression for the absorbed dose can be re-written in the 
form 


D(xsy,Z) = Tf Glx,ys0) F(p,z)do . (7) 
O 


We note that the absorbed dose is an integral over two factors: one of them, 

G, depends only on the beam geometry relative to the point (x,y,z); the other, 
F, 1s the elementary energy-deposition function available from our Monte Carlo 
calculations. Of particular importance is the uniform beam of finite cross 
section, characterized by a beam intensity distribution Q that is constant over 
a certain field size, and zero elsewhere. This will be discussed in Section 5. 


4. BROAD-BEAM DEPTH-DOSE DISTRIBUTIONS 


We now consider the case of a uniform broad beam of unlimited lateral extent, 
i.e., a beam with intensity distribution Q(p,v) = constant. Depth-dose distri- 
butions from such beams (for beam energies between 1 MeV to 60 MeV) are shown in 
Fig. 4 and are given numerically in Table 2. These plots were obtained by drawing 
Smooth curves through Monte-Carlo-generated histograms similar to the histogram 
for a 20-MeV beam shown in the upper right-hand corner of Fig. 1. The depth-dose 
curves are given in scaled form, as plots of the dimensionless quantity (1ro/To)® 
vs. the scaled depth z/ro, with 9M again the energy deposited per unit 

depth. The curves extend beyond the depth z= ro because of the combined 
effect of energy-loss straggling and the transport of energy by secondary 
bremsstrahlung photons. The higher the beam energy, the more noticeable is the 
bremsstrahlung tail of the depth dose curves, which is due to radiation produced 
by electrons within the phantom. 


It is customary to obtain the practical range by fitting a straight line to a 
portion of the depth-dose curve beyond the peak. As illustrated in Fig. 5, one 
can extend the straight line to the depth where the extrapolated dose is zero, 
obtaining the practical range ry; or one can extend the line to the depth 

where it meets the bremsstrahlung tail, obtaining the practical range rx. We 
have analyzed the depth-dose curves in Fig. 4 in this manner, obtaining first 

the ratios r./ro and r5/Vo at several energies, and then the practical ranges 
Ty and rp hith the use of appropriate Yro-values. Linear least squares fits 
were made of the ratios rp/To and rx/To as functions of 1/To;, which lead to 


the following relations (with rp and> r5 in cm of H20, and To in MeV): 


ry = 0), 501. To, =: O10 
-using "old" range values Yo from Table 1 (8) 

mi = 0.496 To - 0.104 

ro =O 505 To. =O 106 
using "new" range values Yo from Table 1 (9) 

a = 0.500 To - 0.101 


It is estimated that the practical ranges have an overall uncertainty of ~ 2%. 


In a further analysis of the depth-dose curves in Fig. 4, we have determined the 
depths Zsq at which the curves have fallen to 50% of their peak value. The 
relationship between Zso and the beam energy To is almost, but not quite, 
linear. In Fig. 6 we give the ratio To/Zso9 as a function of Zsa, and the 
reciprocal ratio Zs9/To as a function of To, from which one can determine 

To as a function of Zs 9 and vice versa. Two sets of curves are given in 

Fig. 6, one (solid curves) obtained with the "new" ranges, and the other (dashed 
curves) with the old ranges from Table 1. 


5. DOSE DISTRIBUTIONS FROM BEAMS WITH CIRCULAR CROSS SECTIONS 


The geometry factor G given by Eq (6) can be evaluated numerically when the 
beam intensity Q is variable, and often analytically when Q has a constant 
value Q. over a prescribed "field size" and vanishes elsewhere. A particularly 
Simple cise is that of a circular field of radius R, for which we shall give 
some illustrative results. Let D(p,z) denote the absorbed dose at a depth z 
in the phantom and at a radial distance p from the axis of the incident beam. 
The expression* for the absorbed dose, obtained by summing over a distribution 

of point-monodirectional beams, is 


D(p.z) = Ty f Go" Elona) Ro 54) ss (10) 
(@) 


By some simple geometric considerations it can be shown that G(p,0') is given 
by the following equations: 


Lf tone O Oe ol = R G =. 
(11a) 

() > 1 G=0 

iner< 0 <-R OS ol S Reso Ges @. 
R-p<p'<R+p G= (Q/m) cos’ '[(p'@+p*-R°)/20'p] (116) 

Rot p =o: G=0 

If p=R 0<p' <2R G = (Q./m) cos” '(p'/2R) 

(ile) 

2? <& oy! G=0 

ip (oy Bel OS oO So oR G=0 
p-R<p'<p+R G= (Q./m) cos [(p'“tp-R“)/2p'p] (11d) 

ae IRS yo)! G=0 


Figure 7 shows central-axis depth-dose curves D(0,z) calculated from Eqs (10) 
and (lla) for beams with circular cross sections. A series of depth-dose curves 
are plotted for a set of increasing beam radii (field sizes), illustrating the 
transition from narrow-beam to broad-beam conditions. The dependence of the 
Shapes of these curves on the field size can be understood in terms of the 
competing effects of electron multiple-scattering into and out of the central- 
axis region. Whereas the outscattering is independent of the field size, the 


“This expression is equivalent to Eq (7); however, x and y are replaced by 
Os aNd oy) 1s replaced by p'. 


inscattering increases with field size, which results in a larger buildup and 
peak and a slower decline beyond the peak for the broad-beam depth-dose curves. 


Figure 8 shows the radial distribution of absorbed dose at various depths, for a 
20-MeV beam with a radius of 1.25 cm, obtained by application of Eqs (10) and 
(1la-1ld). The radial distributions actually are averaged over small depth 
intervals. At shallow depths, the radial distributions are almost constant over 
the irradiated field, and then drop off very rapidly. With increasing depth the 
drop-off becomes less sharp, due to the effect of multiple scattering. The 
radial distributions can be used to generate conventional isodose contours. 
Figure 9 shows isodose patterns derived in this manner for 10- and 20-MeV beams 
with radii of 2.5 and 5 cm. 


6. CONCLUDING REMARKS 


The results in this paper are only illustrative. The body of information con- 
tained in the calculated absorbed-dose distributions at 14 energies from 1 to 

60 MeV remains to be exploited more fully. A very simple and straightforward 
extension would be to compute depth-dose curves and isodose curves for fields of 
various sizes and shapes, including beams with distributed energy spectra. A 
more difficult extension would be to discard the assumption of perpendicular 
incidence and to take into account the angular dispersion of the incident beam. 
Finally, it remains to be investigated to what extent the superposition method 
can be retained when treating inhomogeneous phantoms. 
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TABLE | 


CSDA Range ie in Water 


1 | "0 (/cn’) Percent 
(MeV) o1d@ New? Rae) ea ae 
60 22.65 Zen e 053 
50 19E72 19.79 0.4 
40 16.53 16.62 0.5 
30 13205 13.14 Oey 
20 9.209 9.295 0.9 
hs TENG 7.200 eal 
10 4.902 4.963 le 

8 3.969 4.021 is 

6 3.006 3.048 1.4 

5 (Apso) 4 2.547 1.4 

4 2.008 2.037 1.4 

3 1.493 1.514 1.4 

2 0.9654 0.9802 ee 

| 0.4305 0.4377 Sd) 


Mean excitation energy 65.1 eV; density effect correction 
from Sternheimer and Peierls.° 


Mean excitation energy 75 eV; density effect correction 
from Ashley.” 


TABLE 2 


Scaled distribution of deposited energy per unit depth, for broad parallel 
monoenergetic beams. To obtain the absorbed dose in units of grays (J/kg), 
the entries in this table should be multiplied by 1.602 x 10 °° N To/ro; 
where N is the number of electrons incident per cm*, To is the beam 
energy in MeV, and Yo is the electron range in g/cm’. 
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FIGURE CAPTIONS 
Distribution of energy deposited by a 20-MeY point-monodirectional beam. 


Curves of scaled radial distance ys. depth, along which the normalized 
cumulative radial distribution of deposited energy is constant. 
Fraction of energy deposited at radial distances greater than p/ro 

a. 10-MeV point-monodirectional beam 

b. 20-MeY point-monodirectional beam 


Scaled broad-beam depth-dose curves. 
Definition of practical range. 


Relation between incident beam energy To and depth Zso at which the 
broad-beam depth-dose curve has fallen to 50% of its peak value. 


Central-axis depth-dose curves from beams with circular cross section 
(radius R). 

a. Beam energy 10 MeV 

b. Beam energy 20 Mey. 

Radial distribution of absorbed dose from a 20-MeV circular beam with 


radius R= 1.25 cm. Curves shown are averages over indicated depth 
intervals. 


Isodose curves in a water phantom irradiated with 10- and 20-MeV 
electron beams with circular field sizes. 
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